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Single electron control in n-type semiconductor quantum dots using non-Abelian 
holonomies generated by spin orbit coupling 
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We propose that n-type semiconductor quantum dots with the Rashba and Dresselhaus spin 
orbit interactions may be used for single electron manipulation through adiabatic transformations 
between degenerate states. All the energy levels are discrete in quantum dots and possess a double 
degeneracy due to time reversal symmetry in the presence of the Rashba and/or Dresselhaus spin 
orbit coupling terms. We find that the presence of double degeneracy does not necessarily give rise 
to a finite non-Abelian (matrix) Berry phase. We show that a distorted two-dimensional harmonic 
potential may give rise to non-Abelian Berry phases. The presence of the non-Abelian Berry phase 
may be tested experimentally by measuring the optical dipole transitions. 
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I. INTRODUCTION 

Single electron control in semiconductor quantum dots 
would be valuable for spintronics, quantum information, 
and spin qubits 1]. Adiabatic time evolution of degen- 
erate eigenstates of a quantum system provides a means 
for controlling individual quantum states[2j, |3|, ■ They 
exhibit non-Abelian gauge structures and often give fi- 
nite non-Abelian Berry phases (they are also called ma- 
trix Berry phase or holonomic phase). These phases 
often depend on the geometry of the path traversed 
in the parameter space of the Hamiltonian. Nuclei 0, 
superconducting nanocircuits0, optical systems^, and 
atomic systemsjg, 0, have such degenerate quantum 
states. It has been shown that universal quantum com- 
putation is possible by means of non-Abelian unitary 
operations^! 0]. Manipulation is expected to be sta- 
ble since symmetries of the Hamiltonian that give rise 
to degeneracy are not broken during the adiabatic trans- 
formations. However, the degree of its stability is under 
investigational . 

Recently several interesting possibilities for electronic 
manipulation in semiconductors have been proposed. 
Spin manipulation of quasi-two-dimensional electrons by 
time-dependent gate voltage is possible through the Dres- 
selhaus and Rashba spin-orbit coupling mechanisms |14| . 
A geometric spin manipulation technique based on ac- 
ceptor states in p-type semiconductors with spin-orbit 
coupling was proposed [15|. Spin-orbit coupling and a 
revolving external electric field may generate holonomic 
qubit operations in CdSe nanocrystals 0|. Holonomic 
quantum computation using excitons in semiconductor 
nanostructures has been also proposed|l7j. 

A matrix Berry phase is experimentally interesting 
because it represents mixing between degenerate lev- 
els. Let us explain briefly the basic ideas using a sim- 



ple system possessing Kramers' double degeneracy ||. 
The Hamiltonian depends on some external parameters 
X p . If the system is in a superposition state \ip'(0)) = 
ci (0)|^i (0)) +C2(0)|^2(0)) at time t=0 an adiabatic evo- 
lution of the parameters A p can transform this state into 
another state \ip'(t)} = a(t)\ipi(t)) + c 2 {t)\vb 2 {t)) after 
some time t. Here the orthonormal basis states \ipi(t)) 
are the instantaneous eigenstates of the Hamiltonian. For 
a cyclic change with the period T, represented by a closed 
contour C in the parameter space, the states \ipi(T)) and 
\ip 2 (T)) return to the initial states IV'i(O)) and | V2 (0)} , 
but the coefficients ci(T) and c 2 (T) may not return the 
initial values. In such a case a 2 x 2 matrix Berry phase 
(non-Abelian Berry phase) <&c is generated 



ci(T) 

C2(T) 



ci(0) 
c 2 (0) 



(1) 



This non-Abelian geometric phase (holonomy) is con- 
nected to non-Abelian gauge potentials, as we discuss 
below. During the adiabatic cycle the degenerate energy 
E(t) varies with time. But hereafter we will set E(t) = 
since it can be easily restored when needed. The expan- 
sion coefficients satisfy the time-dependent Schrodinger 
equation 



ihcj 



13^3 



i= 1,2. 



(2) 



The matrix elements Aij are given by A 



h^2 P (Ap)i,j-^f, where the sum over p in Ai 



dt 

to be the sum over A r , 



A p are 2 



is meant 
2 matrices and are 



called the non-Abelian vector potentials. They are given 
by 



{A P )i,i=m^-) 



(3) 
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The formal solution of the time-dependent Schrodinger 
equation gives the matrix Berry phase $c* — 
p e ~ § c £p A p d\ p ^ w j iere p represents a path ordering. 
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FIG. 1: (a) An electric field along the z-axis quantizes the 
electronic motion in a triangular potential along the axis. In 
our work the triangular potential is sufficiently strong and 
only the lowest energy subband is included. The structural 
inversion asymmetry in V(z) leads to the Rashba interac- 
tion. An adiabatic change can be induced by changing the 
magnitude of the electric field, (b) The two-dimensional elec- 
tronic motion is quantized in a distorted harmonic potentials 
U(x,y) = ^m*L0 2 x 2 + ^m*uj 2 y 2 + e'y. An adiabatic change 
can be induced by changing the magnitude of the distortion 
potential, (c) In the presence of the Rashba and/or Dressel- 
haus spin orbit coupling terms each discrete eigenstate of a 
semiconductor quantum dot has a double degeneracy due to 
time reversal symmetry in the absence of a magnetic field. 
The three lowest energy levels are shown schematically (Each 
pair is doubly degenerate), (d) A schematic drawing of a 
cyclic adiabatic path. The parameters En and E p depend on 
the magnitudes of the electric field and the deviation from the 
two-dimensional harmonic potentials. 



Under a unitary transformation \^)[) — ^2jU^\ipj) the 
non-Abelian gauge structure emerges 

A' k = UA k U^+iU^-. (4) 

Here we propose that quantum dots in n-type semicon- 
ductors with spin-orbit interaction can have matrix Berry 
phases. All the discrete energy levels of quantum dots 
possess a double degeneracy because the Rashba and/or 
Dresselhaus spin orbit coupling terms have time reversal 
symmetry. In II- VI semiconductors the Rashba term is 
expected to be larger than the the Dresselhaus coupling. 
In III-V semiconductors, such as GaAs, the opposite is 
true^^- I n our work both the Rashba and Dresselhaus 
terms are included. The adiabatic transformation can 
be performed electrically by changing the confinement 
potentials of the quantum dot, see Fig^d). We find 
that, although the spin orbit terms formally give rise to 
a non-Abelian structure for the matrix vector potentials, 



double degeneracy does not necessarily lead to finite non- 
Abelian Berry phases. We show that when the circular 
symmetry of two-dimensional harmonic potentials is bro- 
ken matrix Berry phases can be produced (such a distor- 
tion is shown in Fig^b)). We propose that the presence 
of a matrix Berry phase may be detected by measuring 
the the optical spectrum. 

In Sec. II we describe the Hamiltonian of the system 
in detail, and in Sec. Ill we discuss a model calculation 
of the matrix Berry phase. A possible experimental de- 
tection of the the matrix Berry phase is suggested in Sec. 
IV. Conclusions are given in Sec. V. 



II. HAMILTONIAN 

An electric field E is applied along the z-axis and elec- 
trons are confined in a triangular potential V(z), see 
Fig^a). When the the width of the quantum well along 
the z-axis is sufficiently small we may include only the 
lowest subband state along the z-axis. We denote this 
wavefunction by f(z), see Fig^a). From the expec- 
tation value (f(z)\h%\f(z)) = 0.8(2m*eE/h 2 ) 2/3 we es- 
timate the characteristic length scale along the z-axis: 
R z = l/^/0.8(2m*eE/h 2 ) 2 / 3 . The Hamiltonian in the 
absence of the spin orbit coupling is Hk = — 6 2 J, + 
U(x,y) + V(z). We take the two-dimensional poten- 
tial to be U(x,y) = \m*Lo 2 x x 2 + \m*<jj 2 y 2 + V p (x,y), 
see FigQIb). The strengths of the harmonic poten- 
tials are denoted by u> x and u> y . They may be varied 
by changing gate potentials of the quantum dot system. 
The characteristic lengths scales along x- and y-axis are 

R x y = '■ The potential V p (x,y) — e'y represents 

a distortion of the harmonic potentials and its strength 
e' may be varied electrically. 

In a periodic crystal potential of a semiconductor the 
spin orbit interaction has two contributions. The Rashba 
spin orbit term is 

Hr = CR, (a x ky - (Jykx) ■ (5) 

Here a XjV are Pauli spin matrices and k xy are momen- 
tum operators (k x = A-|- and similarly with k y .). The 
constant cr depends on the external electric field E ap- 
plied along the z-axis. The Dresselhaus spin orbit term 
is 

H D = c D ((a x k x (k 2 y - k 2 z )) + (a y k y (k 2 z - k 2 x ))) . (6) 

There is another term of the form <j z (k z ) [k 2 — k 2 ) in the 
Dresselhaus spin orbit term but it vanishes since the ex- 
pectation value (k z ) — (f(z)\k z \f(z)) = for the first 
subband along z-axis. The constant cd represents break- 
ing of inversion symmetry by the crystal in zinc blende 
structures. 

The total Hamiltonian of an electron in a semiconduc- 
tor quantum dot is H = Hk + Hr + Hr> ■ An eigenstate 
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of the Hamiltonian H may be expanded as a linear com- 
bination of the eigenstates of Hx 



E 



Cmn^mn T) + > y Cmnl\mn I). (7) 



The expansion coefficients satisfy a matrix equation 
J2 m 'n'a'(' mna \ H \ rn ' n ' al ) c m'n'a' = Ec mna . In the basis 
states \mna) the quantum number m(n) and a label the 
harmonic oscillator levels along the x-axis (y-axis) and 
the component of electron spin. The subband wavefunc- 
tion f(z) is suppressed in the notation \mncr). 

In the absence of the Zeeman term the total Hamilto- 
nian is invariant under time reversal symmetry: k — > — k 
and a — *■ —a. The time reversal operator is K — —ia y C, 
where the operator C stands for complex conjugation. 
The time reversed state of is 



E 



C*nnl\ mn T) 



mn 



mn[). (8) 



Note that K 2 \tfj) = —\tp). These two states are degen- 
erate and orthonormal. We have suppressed the Bloch 
wavefunction of the conduction band in applying the time 
reversal operator since it is unaffected by the operator 
K. Our wavefunctions are all effective mass wavcfunc- 
tions and only the conduction band Bloch wavefunction 
at k — is relevant. 

An adiabatic change implemented by changing the en- 
ergy parameters Er and E p that characterize the Rashba 
term and the distortion potential V p {x,y): 



Eft = 



CR 

V2Ry 



and E v = eR v 



(9) 



where e = e' /V2. The first parameter Er depends on 
the electric field through the constant cr. The typi- 
cal value of the energy scale associated with the Rashba 
constant depends on the electric field applied along the 
z-axis and the semiconductor material: it is of order 
Er = cr/R ~ 0.01 — lOmeV, where the length scale 
R ~ lOOA is the lateral dimension of the quantum 
dot. The second adiabatic parameter E p represents the 
strength of the distortion potential e'y: the expectation 
value of the distortion potential is E p = (Q\e'y\i). Its 
magnitude is of order 1 — lOmeV, depending on the elec- 
tric field applied along the y-axis and the width of the 
triangular potential V(z). 



III. MODEL CALCULATION OF 
NON-ABELIAN BERRY PHASE 

A. A truncated Hamiltonian matrix 

We work out a model that can be solved analytically. 
This model is simple but much can be learned from it. 
Let us take lo x = 2tv y (Other values of lo x can also be 
chosen). Then the lowest eigenenergy state of Hk is 



\mn) — 1 00) with the energy Eq = '\hw y and the next 
lowest eigenenergy state is |01) with the energy Ei = 
|_E - The typical value of the energy spacing between 
the quantum dot levels, Eo, is of order 1 — lOmeV. Let 
us write down the eigenstates of the total Hamiltonian 
as a linear combination of four basis states made out of 
these states and spin degree of freedom: 

|V) = co,o, t |0,0,T) + co,i,t|0,1,T) 

+ co,o,i|0,0,l) + c , u |0,l,|)- (10) 

The 4x4 truncated Hamiltonian matrix is 



Eo 


E v 





-iE R - E D 


E p 


Ex 


iE R + E D 








-iE R + E D 


Eo 


Ep 


iE R - E D 





Ep 


Ei 



(11) 



or 



±(E + E 1 )I+±(E -E 1 ] 1 " 



a z 



-E n 



-E D 



cr x 
a x 

- 



a v 
a y 



(12) 



The matrix elements of the distortion potential are 
(m|y|n.) = J ' 2n [i UJy (Vn + I5 m . n+ i + \fn5 m ,n-i)- Note 
that the distortion potential couples even and odd parity 
states of the one-dimensional harmonic potential. Other 
functional form of V p (x,y) may be also used to gen- 
erate the non-Abelian Berry phase as long as it cou- 
ples even and odd parity states. The first term in 
the Dresselhaus spin orbit term, Eq. ©, is zero since 
(0|fc 2 |0) = 0. From the second term of the Dresselhaus 
term we get -ic D <0|^|l)(</(*)|^|/(*)) - <0|fc*|0)) = 
—Ed, where the constant Ed = I{Er) — g(E ) with 
f(E R ) = c D /R y R 2 z and g(E ) = c D /R y Rl- Here we 
have used the momentum matrix elements (m\hk x y \n) = 



-(\Ai + l^m,n+i - y/nSm.n-i)- The magnitude 



of the energy scale associated with the Dresselhaus term 
is of order Ed = cd/R 3 , and it can be larger or smaller 
than the Rashba term, depending on the material |T^. 
The function /(Er) depends on Er because the electric 
field E enters through R z . There is no simple way to 
calculate cr because it depends both on the electric field 
inside the semiconductor heterostructure and on the de- 
tailed boundary conditions at the interface. For simplic- 
ity we take /(Er) — ciEr and g{Eo) = bE$, where a and 
b are numerical constants. More complicated functions, 
for example, J(Er) — E R /E and g(E ) = E , could 
be used, but our calculation indicates that the essential 
physics does not change. 
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B. Finite matrix Berry phase 



Diagonalization 
matrix gives the 



El 



9E% 



9Ej 



of this 4 
eigenenergies 

Let 



9E R ). 



x 4 Hamiltonian 
A = ±(4£ ± 

us choose the fol- 



lowing state 



3E n 



En 



El 



9(E% + El 



V 



D 

3(E D - iE R 




El) 



(13) 



and its time reversal state 
/ 



1 

~N* 



3(E D + iE R ) 

3E„ 



\E - JEl + 9(E 2 D +E$+E%) J 



(14) 



as the basis states in the lowest energy degenerate Hilbert 
space. 



of 



the matrix 

|2, 



An 



we use 



In the evaluation 
2J^W*%P = £;f\Mr)\ 2 dr = 0, where k 
mna. Note that </>fc(r) is the two-dimensional harmonic 
oscillator wavefunction and that it is a real function. 
One can show that (A p )ij = i J^ fc o^^^S where the pair 

of degenerate states are = J2k a k\k) an d \^j) = 
"Y^uPk\k)- The orthonormalization (i/>i\ipj) — §ij gives 
that the diagonal matrix elements (A p )i^ are real and 
that the off-diagonal elements satisfy (A p )ij = (Ap)*^. 
The adiabatic parameters are Ai = Er, A2 = E p . We 
calculate the matrix Berry phase for Ed = Er — En. 
The matrix vector potentials with respect to \ip) and \ip) 
have the following structures 



A Er = 



ci ai + ib\ 
a\ — ib\ —C\ 



(15) 



and 



A» = 



a-z — ib 2 
a 2 + ib 2 

a 2 a x + b 2 a y . 



(16) 



The functions a,, bi, and Ci depend on Er and E p and 
are real. Under an adiabatic time evolution a state 
in this degenerate Hilbert space changes as \ip'(t)) = 
cx(t)\i/j(t)) + c 2 (t)\ip(t)). Suppose initially ci(0) = 1 and 
c 2 (0) = 0, i.e., |^'(0)) = The cyclic adiabatic path 
with the period T = 2it/lu is shown in Fig^d) and is 
given by 



( E R (t),E p (t)) = 

( E R . C + AE R cos(wt), E Pj , 



• AEpSin(wi)). (17) 



The frequency u> can be taken to be a fraction of En. We 
solve Eq.J2J numerically and find that c 2 (T) is non-zero. 
For the parameters Er, c = 2E , E p _ c = En, AE R = 
1.9E , AE p = 0.9E , and us = En/10 we find c x (T) = 
0.8884 - i0.0897 and c 2 (T) = -0.4429 + i0.0874. 



C. Absence of matrix Berry phase 

The matrix Berry phase is absent when the distortion 
potential V p (x,y) is zero. When only Rashba and/or 
Dresselhaus terms are present in Ea. pi|l the adiabatic 
transformation can be performed by varying the param- 
eters Ai = Er and A2 = En. In this case an orthonormal 
basis set exists in the degenerate Hilbert space such that 
the matrix vector potentials are diagonal and the non- 
Abelian Berry phase is zero. This degenerate basis set, 
\ipi) and \ip2), has the property that for each fc either 

= or e£' = 0, where \ipi) = J2k=i c k\k) ( see ' f° r 
example, Eas. (|18|) and $19$ ). The off-diagonal matrix 



elements of the vector potential (A p )i <2 — ii^iljx) 



■ V 4 (l)*dc 
1 L~ik=\ c k 



because either c 



(i) 



or c 



(2) 



0. 



(2) 

d\p k k 

Therefore the matrix vector potentials are diagonal. In 
this case the matrix Berry phase will be absent. Let us 
construct explicitly \ipi) and \ip 2 ) when only the Rashba 
term is present, i.e., when Ed = E p = 0. They are given 
by 



|</>i 



V2\/El 



VE R 



E Q ^/Ej+9E 

( i(E Q + ^JE> + 9E R ) \ 


3E R 



J 



(18) 



and 



W = 



Eo-VW+VEr- 

\ 



i(E - y/E?£ + 9E R ) 
3E R 




(19) 



With these new eigenstates it is possible to show that the 
non-Abelian gauge potentials are not only diagonal but 
also zero: Ae — and Ae r = 0. 



D. Non-Abelian gauge structure 

Let us also test the non-Abelian gauge structure given 
by Eq. J3J. We make a unitary transformation from 
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(IV'i), 1^2)), given in Eqs.^SJ) and i|19|l . to a pair of time 
reversed degenerate eigenstates 



( 



W) = 



and 



2y/E$+9E[ 



Er 
Er 



\ 



\i[E + ^E* + QEl 



(20) 



9^1 




(21) 



Abelian gauge potentials are 



A' 



1 



El + 9El 



3# /2 
3E /2 



3En 



2{E% 



9E%) 



and 



A' 



1 



EI+9EI 
ZE R 





-ZEr/2 



-3E r /2 




2{El 



9E\) 



(23) 



Since the old vector potentials are Ae r 



it follows from Eq. (gj that A' k 



iU 



and A Eo = 
. We have ex- 



plicitly verified that this relation holds by computing the 
unitary matrix U. We have also verified independently 
the absence of the matrix Berry phase by solving the 
time-dependent Schrodinger Eq.(J5J with A' Er and A' Eq . 
This provides a check on our numerical method of solving 
the time-dependent Schrodinger equation. 



IV. DETECTION OF MATRIX BERRY PHASE 

After an adiabatic cycle the electron acquires a ma- 
trix Berry phase. The presence of such a phase may 
be detected by measuring the strength of dipole optical 
transitions before and after the cycle. In this section 
we calculate the optical strengths using the truncated 
Hamiltonian. This calculation is not quantitatively ac- 
curate but it will give us an estimate of the magnitude 
of the effect. 

First we need to prepare physically some particular 
pair of degenerate states. In the presence of a magnetic 
field along the z-axis any degenerate pair of states at zero 
magnetic field will be split into two states. We define the 



lowest energy pair of degenerate eigenstates in the zero 
magnetic field limit of B z as \ipi) = ]ira.B x -*o |"0i(-B z )) 
and \ip 2 ) = lim Bz ^ \ip 2 (B z )) , where \^i{B z )) and 
\4 ) 2{B Z )) are the split lowest and next lowest energy 
states in a finite magnetic field, see Fig.(|5Ji). In a sim- 
ilar way we define degenerate eigenstates in the zero 
magnetic field limit of B x : |0i) = lim^^o \<Pi{B x )) 
and 1^2) = lims x ^o \4>2(B X )), see Fig.(|2)). The states 
|'0i,2(-Bz)) and \4>i,2 (B x )) are calculated from the Hamil- 
tonian by replacing k with k+^A and adding the Zeeman 
term (A is the vector potential and e > 0) . 




l¥ 2 (B z )> 
l¥i(B z )> 




I 4> 2 (B X )> 



non- 


B 

(a) * z 

E 


B 

(b) * x 


(22) 






•Kx^^^y^ photon 




l4' 1 (0)>,l^ 1 (T)> 




l¥,(B 7 )> 









(c) Bz 



FIG. 2: (a)|^i) and \1jj2) are generated in the zero limit of 
B z . (b)|<^>i) and \(j>2) are generated in the zero limit of B x . 
(c)Dipole transition from the first to third lowest energy states 
in a finite B z . The initial state in this transition is |</)i(0)) 
or \<j)\ (T)). These states represent, respectively, the states 
before and after the adiabatic cycle. The incident photon is 
polarized along the y-axis. 



Suppose that the electron is in the lowest energy state. 
In order to measure the matrix Berry phase we perform 
the following set of procedures: 

1. We apply a magnetic field along the x-axis and 
take the zero field limit. The resulting state |0i) 
(Fig.JSJs)) can be written as a linear combination of 

and |^a) (Fig. 00): |0i)=ci(O)|Vi)+c 2 (O)|^ 2 ) 
with (| Cl (0)| 2 ,|c 2 (0)| 2 H(l/2,l/2). 

2. We apply adiabatically a magnetic field along the 
z-axis. We increase it to a value B z . Then the prob- 
ability for the electron to be in the state |V>i(-B 2 )) 
is given by |ci(0)| 2 . 

3. We then measure the intensity of the dipole tran- 
sition from the lowest energy state to the third 
lowest energy state, see Fig.(|3;). This intensity 
is proportional to |ci(0)| 2 , i.e, proportional to the 
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probability that \ip\(B z )) is occupied. The lowest 
and third lowest energy states can be written as 
a linear combination of the basis states \mna) : 
both states can be written in the form cq o ||0, 0, | 
> + c , 1 , T |0, 1, T)+co,o,||0,0, O+co^lO, 1, i)'. When 
the photon is polarized along the y-axis only the ba- 
sis states with different parities are coupled in the 
dipole approximation, for example, (00cr| fc y |01cj) = 
— (01cT|A; y |00o") is non-zero. 

Again we assume that the electron is in the lowest energy 
state. We perform the second set of procedures: 

1. We apply a magnetic field along the x-axis and take 
the zero field limit. 

2. We then perform an adiabatic cycle following a 
closed loop in the parameter space of (Er,E p ), 
given by Ea. 1)17(1 . After the adiabatic cycle the 
electron will be in the state \<pi(T)) = ci(T)|^i) + 
c 2 (T)|V 2 ) with (| Cl (r)| 2 ,|c 2 (T)| 2 ) = (0.559,0.441) 
(the parameters are Er <c = Eq, E p>c = 3Eq, 
AE R = 0.8-Eo, AE p = 2.5E , and uj = 6.2E ). 

3. We apply adiabatically a magnetic field along the z- 
axis, see Fig. 10;) • Then the probability for the elec- 
tron to be in the first lowest energy state, \tpi(B z )}, 
is |ci(T)| 2 while the probability that the electron 
will be in the second lowest energy state, |?/> 2 (-B z )), 
is |c 2 (T)| 2 . Here the value of B z is the same as that 
of step 2 in the first set of procedures. 

4. We then measure the intensity of the dipole transi- 
tion from the lowest energy state to the third low- 
est energy state. This intensity is proportional to 
|ci(T)| 2 , i.e, proportional to the probability that 
\ifJi{B z )} is occupied. 

These two sets of measurements are repeated many 
times. Then the ratio i^q)^ = 1.12 gives the intensity 
ratio of the dipole transitions in the first and second sets 
of procedures (the dipole matrix elements cancel each 
other). It is the direct measure of the matrix Berry phase. 
The analysis of the intensity of optical transitions in zero 
magnetic field is complicated due to the presence of the 



degeneracy in the final states of the transition. In a finite 
magnetic field this degeneracy is lifted. 



V. DISCUSSIONS 

Each discrete eigenstate of a semiconductor quantum 
dot with the Rashba and/or Dresselhaus spin orbit cou- 
pling terms possesses a double degeneracy due to time re- 
versal symmetry. We have investigated the matrix Berry 
phase of such a quantum dot in a simple truncated model 
Hamiltonian that can be solved analytically. We have 
found that the double degeneracy does not necessarily 
lead to a finite non-Abelian Berry phase. The addition 
of a parity breaking distortion potential to the Hamilto- 
nian when both the Rashba and Dresselhaus spin orbit 
coupling terms are present gives rise to a finite matrix 
Berry phase. We have proposed that this phase may be 
detected in the dipole transitions between the ground and 
first excited states in a magnetic field. 

For an accurate modeling of possible experiments the 
number of basis vectors in the truncated Hamiltonian 
matrix must be chosen sufficiently large. This number 
will be determined by the ratios Eq/E^d^. Calculation 
of the matrix Berry phase in such a case requires a heavy 
numerical computation. Accurate experimental determi- 
nation of the functional dependence of the Rashba con- 
stant on the electric field would be also valuable in de- 
termining the magnitude of matrix Berry phases. The 
quantum dot studied in this work contains a single elec- 
tron. It may be worthwhile to investigate the effect of 
many body interactions. 
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